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A microscopic model is proposed for the scalar order in PrFe4Pi2 where crystalhne electric field 
(CEF) singlet and triplet states interact with two conduction bands. By combining dynamical mean- 
field theory and the continuous-time quantum Monte Carlo, we obtain an electronic order with staggered 
Kondo and CEF singlets with the total conduction number being unity per site. The ground state becomes 
semimetallic provided that the two conduction bands have different occupation numbers. This model nat- 
urally explains experimentally observed properties in the ordered phase of PrFe4Pi2 such as the scalar 
order parameter, temperature dependence of the resistivity, field-induced staggered moment, and inelastic 
features in neutron scattering. The Kondo effect plays an essential role for ordering, in strong contrast with 
ordinary magnetic orders by the RKKY interaction. 
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Among Pr-based filled skutterudites that have attracted 
much recent attention,'^ PrFe4Pi2 is one of the most inter- 
esting; it shows a clear Kondo-like behavior in the resistivity, 
and undergoes a second-order phase transition at Tc = 6.5K 
into non-magnetic staggered order In spite of the clear cusp 
in the magnetic susceptibility^- at Tc there is no magnetic 
moment in the ordered phase. Instead, the neutron diffraction 
has revealed the field-induced staggered magnetic moment,"^' 
which shows the nonmagnetic character of the staggered or- 
der parameter. Theoretical^'*^ and experimental^^ studies have 
already identified the order parameter as a scalar, but the mi- 
croscopic nature of this mysterious scalar order is yet to be 
clarified. 



The inelastic neutron scattering experiments^"'"' on 
PrFe4Pi2 have probed the broad quasi-elastic response above 
the transition temperature, which is characteristic of the 
Kondo effect. Below Tc, on the contrary, sharp inelastic ex- 
citations have been observed in addition to remaining quasi- 
elastic components. The Uneshape of the spectrum does not 
depend significantly on wave vector,"* but the intensity of 
the scattering is stronger at the center of Brillouin zone than 
at the boundary.'"' The observed Kondo behaviors indicate 
a tendency toward itinerant /-electron state, while the in- 
elastic feature may be interpreted as crystalline electric field 
(CEF) excitations. Hence, PrFe4Pi2 can be regarded as a sys- 
tem located on the boundary between itinerant and local- 
ized characters of electrons. Accordingly the ordered phase 
in PrFe4Pi2 may not be understandable in terms of the con- 
ventional RKKY interaction. 

In the present paper, we propose microscopic mechanism 
for the scalar order in PrFe4Pi2. In our previous work,'^' we 
found a new electronic order with staggered Kondo and CEF 
singlets in a generalized Kondo lattice. However, the previ- 
ous paper obtained an insulating ground state, and did not 
include comparison with experimental results in PrFe4Pi2. 
In this paper, we demonstrate that the staggered Kondo-CEF 
singlet order is a plausible scenario for PrFe4Pi2 to explain 
its characteristic behaviors qualitatively. For example, we ob- 
tain semimetallic ordered state by taking a modified model 
for conduction electrons. As a result, temperature dependence 
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of the resistivity in PrFe4Pi2 is reasonably reproduced in the 
whole temperature range. Furthermore, sharp inelastic peaks 
observed only below the transition temperature is naturally 
explained as local excitations in CEF sites generated by the 
ordering. 

We take the singlet-triplet CEF states for low-lying 
states in PrFe4Pi2 as have been suggested in ref.l3. For con- 
duction bands, we adopt different orbitals 7 — 1,2 since 
PrFe4Pi2 have two bands composed of p- and d-electrons.''*' 
We take the following model: 



n 



(1) 



where Sc-yi denotes the spin of conduction electrons at site i. 
The local singlet-triplet states are represented by two pseudo 
spins Si and S2, and the CEF splitting A is simulated by 
the coupling between them as shown in the third term. The 
second term shows the Kondo exchange interaction between 
localized and conduction spins. 

We use the dynamical mean-field theory (DMFT) formu- 
lated in the two-sublattice system in order to discuss the 
staggered order. In this theory, the lattice is mapped to the 
effective impurity problem with full on-site correlations.'''' 
The continuous-time quantum Monte Carlo method (CT- 
QMC)'*"^' is applied as the impurity solver. For analytic con- 
tinuation from imaginary axis onto real one, we employ the 
Pade approximation. 

For the spectrum of each conduction band, we take Efci = 
£k + Sfj,/2 and £^2 = £k — Sii/2 with Sfi controlling the 
asymmetry. Here gives the semi-circular density of states 
(DOS) defined by 

2 



Po(£) 
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The band width is given by 2D with = 1 as a unit of en- 
ergy. Throughout this paper, the number of conduction holes 
per site is fixed at ric — ni+n2 = 1, which corresponds to the 
3/4 filling of two conduction bands. In this paper, we take the 
parameters such as J = 0.8, A = 0.2 and (5/i = 0.1. With- 
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Fig. 1. (color online) Density of states in (a) disordered phase at T = 0.04 
and (b) ordered phase at T = 0.01. The inset in (a) illustrates the bare 
DOS of conduction bands with different scales. 




Fig. 2. Temperature dependences of the resistivity R both above and below 
the transition temperature Tc ~ 0.031. The inset shows the low tempera- 
ture behavior For comparison, the ferromagnetic case with J = —0.8 is 
also shown. 



out interaction, we obtain ni ~ 0.56 > n2 — 0.44 at zero 
temperature. Note that the ratio ni/n2 depends on temper- 
ature. The DOS without interaction is schematically shown 
in the inset of Fig. 1. These bands simulate the unit num- 
ber of conduction holes in LaFe4Pi2,''^^ and make the start- 
ing point for PrFe4Pi2 with localized picture of / electrons. 
The Kondo temperature for the impurity system is given by 
Tk = Dex.p{-2D/J) ~ 0.082, which is comparable to A. 
The large value of J is taken simply for easier numerical cal- 
culation. Because of strong renormalization effect, however, 
the characteristic temperatures are still much smaller than the 
band width 2D. Hence, the low-energy behaviors are not sen- 
sitive to the choice of the bare parameters. 

Under the present condition, the system undergoes the stag- 
gered Kondo-CEF singlet ordered phase below T — ^ 
0.031 by the second-order transition. In this phase, the con- 
duction holes accumulate more on the Kondo-singlet site than 
on the CEF-singlet site.'^^ Clearly the symmetry of the order 
parameter is a scalar The origin of this ordering can be under- 
stood from the strong coupling limit with respect to J.^""^^^ 
In this Umit, each site forms the Kondo singlet in presence of 
one conduction electron, or CEF singlet with fully-occupied 
conduction state. The second-order perturbation with respect 
to the hopping of conduction electrons gives the effective re- 
pulsive interaction between Kondo singlets, which favors the 
staggered order with Kondo and CEF singlets. This order is 
induced by the Kondo effect, and the RKKY interaction is 
irrelevant here. As is clear from the formation of the CEF sin- 
glet, the present order is characteristic of the systems with 
P configuration. In terms of itinerant and localized charac- 
ters, the staggered Kondo-CEF singlet order can be regarded 
as alternating sites with itinerant and localized / electrons. 
We note that the present order has nothing to do with nest- 
ing property of the half-filled single conduction band, which 
can be relevant in other skutterudites such as PrRu4Pi2.''*' In- 
stead, it is essential that we have two conduction bands each 
of which is 3/4-filled on the average. 

Now we discuss the physical properties of the order and 
the possible relevance to the experiments in PrFe4Pi2. Figure 
1 shows the DOS both in (a) the disordered phase above Tc 
and (b) the ordered phase. For T > Tc as shown in Fig. 1(a), 



the DOS at the Fermi level is finite. The shape of DOS is mod- 
ified much from the bare one shown in the inset because of the 
strong interaction effect. At T = 0.04, the difference of num- 
bers of holes between two bands is given by ni—n2 ~ 0.094. 
Below the transition temperature as shown in Fig. 1(b), on 
the other hand, a pseudo-gap structure arises near the Fermi 
level with a sharp peak at each edge. The peaks above the 
Fermi level come from Kondo-singlet sites, and the ones be- 
low are due to the formation of CEF-singlet sites. In the or- 
dered phase, the numbers of conduction holes come closer to 
each other with ni - ^2 ~ 0.004 at T = 0.01. The DOS at 
the Fermi level is tiny but finite, and hence the ground state 
becomes metallic. 

Next we derive the electrical resistivity R from the current- 
current correlation function. According to the result for the 
Hubbard model in infinite dimensions, ^-^^ we take the follow- 
ing form for the cuiTent operator J for the two-sublattice sys- 
tem: 

^ = X! X! Wfc7(4A7^CfcB7^ + h.c), (3) 
k 70" 

where Ckx-^a is the annihilation operator of conduction elec- 
trons at sublattice A = A, B with spin a and orbital 7. The 
form of J derives from the nearest-neighbor hopping term 
of conduction electrons in the bipartite lattice, and is off- 
diagonal with respect to sublattice index A. The summation 
is taken over the reduced Brillouin zone of the superlat- 
tice. Then we calculate the quantities 

n{iv„,)= / dT(T,^(T) J) exp(iz.™T) (4) 
Jo 

R-^ = A lim Imn(w + \5)luj (5) 

where Vm = irmrT is a bosonic Matsubara frequency and 
A is a constant. In the framework of the DMFT, the vertex 
correction can be neglected,^'*' and we evaluate the simple 
particle-hole bubble in eq.(5). Assuming that the relevant en- 
ergy range is much smaller than the band width, we replace 
the velocity v^-^ by a constant. 

Figure 2 shows the electrical resistivity R. We also show 
the result with J = —0.8 for comparison, where the Kondo 
effect does not operate. At high temperatures, the resistivities 
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with J = 0.8 and —0.8 tend to the common value obtained by 
the Bom approximation. The resistivity for J = 0.8 increases 
with decreasing T, while it is not the case for J = —0.8. 
Hence the increase in the resistivity at T > Tc is indeed due 
to the Kondo effect. Below the transition temperature Tc, on 
the other hand, the resistivity shows sharp increase, and then 
decreases after showing the maximum around T ^ ClOTc. 
The sharp increase reflects the formation of energy gap at the 
Fermi level, while the decrease shows the metallic behavior 
at sufficiently low temperatures. The inset in Fig. 2 magni- 
fies the low temperature regime with T < 0.006, where the 
Fermi liquid behavior i? cx is found. Experimentally, the 
resistivity in PrFe4Pi2 increases more abruptly^- just below 
Tc than shown in Fig. 2. This may be due to the short-range 
fluctuation of Kondo and CEF sites, which causes the disorder 
scattering, but is not included in our two-sublattice DMFT 

We now discuss magnetic response of the system. In filled 
skutterudites with the CEF symmetry T/j, the triplet states 
consist of Unear combinations of two triplets under the cubic 
symmetry: r4 and Fs states. With use of the mixing param- 
eter d between them, we write the first-excited triplet states 
as |t) — ^/T~~(P\V^) + d|r4). The dipole operator is repre- 
sented as M — aiSi+ a2S2 with ai, 02 being a function^*"' 
of d. The Curie constant is given by C = (a^ + a^)/'^- Ac- 
cording to ref. 27, the triplet in PrFe4Pi2 is mainly composed 
of r4. Hence we choose d = 0.8 in this paper, which gives^*' 
ai = 3.29 and 02 = -0.85. 

Let us define partial magnetic susceptibilities by 



Xm' {qAVrn) 



dT{TrM^{q, t)M^ i-q)) exp(ii.„,T) 



(6) 



where denotes the magnetic moment at sublattice A. The 
wave vector q in eq. (6) belongs to the reduced Brillouin 
zone. The uniform and staggered susceptibilities in the origi- 
nal Brillouin zone are then given by 
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sgnASgnA'XM (0,ii^m), (8) 



with sgnA = 1 and sgne = — 1. We also introduce the local 

loc,A 



susceptibility Xm 



for sublattice A, and their average by 

[XMiqAl'm) 
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)+Xm (ll^r, 



(9) 



where N is the number of sites. The susceptibilities (7), 
(8) and (9) are calculated from local two-particle correlation 
functions in the two-sublattice system. '^•^'^^ 

Figure 3 shows the numerical results for static suscep- 
tibilities. The uniform susceptibility xm{0,0) shows the 
clear cusp at Tc, as in antiferromagnetic transition. Note 
that the experimentally observed magnetic susceptibility also 
shows a cusp.^'^** On the contrary, the staggered suscepti- 
bility xm{Q,0) does not show any distinct anomaly at the 
transition temperature. The magnitude of xm{Q, 0) is much 
smaller than the uniform susceptibility. Namely, the staggered 
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Fig. 3. (color online) Temperature dependences of static magnetic re- 
sponse: Curie law, uniform, and staggered susceptibilities. The field- 
induced staggered susceptibility xs is defined in the text. The inset shows 
local susceptibilities for Kondo- and CEF-singlet sites, together with their 
average. 
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Fig. 4. (color online) Local magnetic spectrum at several temperatures. 



scalar order has nothing to do with the antiferromagnetic cor- 
relation. The inset shows the local susceptibility that splits 
into two different values between Kondo- and CEF-singlet 
sites. The CEF-singlet site has larger susceptibility than the 
Kondo-singlet site. 

We can also discuss the field-induced staggered moment 
observed in the neutron diffraction.^' The relevant susceptibil- 
ity is the staggered magnetic response against uniform mag- 
netic field. Namely, we introduce 



Xs 



(10) 



which vanishes in the disordered phase. Here we have used 
the relation xjiP = xf/^- As shown in Fig. 3, we obtain fi- 
nite Xs in the ordered phase, which means that the staggered 
magnetic moment appears under uniform magnetic field. The 
magnitude of the moment is 61% of the uniform moment for 
T — 0.01 under the present condition. 

The temperature dependence of the local magnetic spec- 
trum Imx^/(cj + i5)/ijj is shown in Fig. 4. Above Tc, the 
spectrum shows broad quasi-elastic peak which is character- 
istic of the Kondo effect, even though the CEF splitting is 
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present. In the ordered phase, on the contrary, the sharp in- 
elastic peak at w ~ 0.14 grows with decreasing temperature. 
This peak is due to emergence of the CEF-singlet site, and 
the corresponding energy is smaller than the magnitude of the 
original CEF splitting, owing to the interaction effect. The 
other and wider inelastic peak at w ~ 0.19 originates from 
the Kondo-singlet site where the DOS has the gap structure. 
These results are consistent with neutron scattering results^^ 
in that the broad quasi-elastic peak is observed in the disor- 
dered phase, and sharp inelastic peaks appear in the ordered 
phase. 

The neutron scattering on single crystal' shows that the 
inelastic peak has almost no dispersion in the ordered phase. 
On the other hand, the intensity of the spectrum depends on 
the wave vector; the intensity is the largest at the zone center 
and is the weakest at the zone boundary. Note that the sum 
rule of the inelastic spectrum with fixed wave number gives 
the static susceptibility. As dynamical magnetic response, we 
have derived only the local spectrum in this paper. How- 
ever, the energy-integrated neutron spectrum can be compared 
with the wave- vector dependent static magnetic susceptibility. 
From the results shown in Fig. 3, it is clear that the dependence 
on wave vector of the neutron intensity is consistent with that 
of the static susceptibilities. 

Let us comment on other interesting aspect of PrFe4Pi2 
in the light of the present results. Since the formation of the 
Kondo singlet between conduction and / electrons is essen- 
tial, the order should be fragile against disorder. In fact, slight 
substitution of Pr by La easily destroys the order as observed 
in Pri_;j:Laa;Fe4Pi2.^'' On the other hand, it has been reported 
that antiferromagnetic insulator is realized under high pres- 
sure. ^'^'^'' This insulating behavior implies the nesting prop- 
erties of conduction bands as in PrRu4Pi2.'^^ We ascribe this 
behavior to the change of the conduction bands. Namely, if 
one of the conduction bands disappears under pressure, the 
remaining band is half-filled, and strong nesting property may 
emerge. Then it is not surprising that we have the insulating 
ordered phase with antiferromagnetism. 

Some experiments show the need of including higher CEF 
levels in PrFe4Pi2. Under high magnetic field along (1,1,1) 
direction, new ordered phases and non-Fermi liquid behav- 
ior have been found. The special property with the (1,1,1) 
field cannot be understood in the singlet-triplet CEF levels, 
but the full CEF model reproduces the level crossing only 
along (1, 1, 1) field.^' Some refinement of the model such as 
the quasi-sextet CEF model including higher doublet will im- 
prove the correspondence between theory and experiments. 

In conclusion, we have demonstrated the relevance of a 
variant of the Kondo lattice model to PrFe4Pi2 using the 
DMFT+CT-QMC. Our scenario of the staggered Kondo-CEF 
singlet order can naturally explain many experimental results 
in PrFe4Pi2 such as the scalar order parameter, temperature 
dependence of the resistivity, and sharp inelastic magnetic re- 
sponse observed only in the ordered phase. 
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